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ABSTRACT 

We present an analysis of interstellar absorption along the line of sight to 
the nearby white dwarf star HZ43A. The distance to this star is 68±13 pc, 
and the line of sight extends toward the north Galactic pole. Column densities 
of Ol, N I, and Nil were derived from spectra obtained by the Far Ultraviolet 
Spectroscopic Explorer [FUSE), the column density of Dl was derived from a 
combination of our FUSE spectra and an archival HST GHRS spectrum, and 
the column density of H I was derived from a combination of the GHRS spectrum 
and values derived irom EUVE data obtained from the literature. We find the 
following abundance ratios (with 2 a uncertainties): Dl/Hl = (1.66±0.28) x 10 1 * * * 5 , 
Ol/Hl = (3.63 ± 0 84) x 10 -4 , and Nl/Hl = (3.80 ± 0.74) x 10 -5 . The Nil 
column density was slightly greater than that of N I, indicating that ionization 
corrections are important when deriving nitrogen abundances. Other interstellar 
species detected along the line of sight were C II, Cm, Ovi, Sin, Ar I, Mgn, and 
Fe II; an upper liinil was determined for N III. No elements other than H I were 
detected in the stellar photosphere. 

Subject headings: ISM: abundances— ISM: clouds— stars: individual (HZ 43A) 
stars: white dwarfs 


1. Introduction 

Deuterium is one of the primary products of big bang nucleosynthesis (BBN), and its 

primordial abundance provides a sensitive measure of the baryon density of the universe 

(Schramm & Turner 1998; Buries 2000). However, deuterium is consumed in stars far more 

quickly than it is produced by the P-P cycle or by any other known stellar nucleosynthetic 
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process (Epstein et al. 1976); thus, measurements of the deuterium abundance at the present 
epoch can only provide lower limits to the primordial abundance. Measurements obtained 
at a variety of redshifts will sample the chemical evolution of the universe, and may permit 
extrapolation to a value for the primordial abundance of deuterium. Measurements of the 
abundance of deuterium and various products of stellar nucleosynthesis, such as oxygen, 
within our Galaxy will improve our understanding of chemical evolution as well as provide a 
lower limit to the primordial deuterium abundance. Recent reviews of deuterium abundance 
measurements can be found in Lemoine et al. (1999), Linsky (1998), and Moos et al. (2001). 

Linsky (1998) reviews measurements of the deuterium abundance for material in 
the local interstellar cloud (LIC), which are consistent with a common value for D/H of 
(1.5 ± 0.1) x 10 5 . The LIC and other nearby clouds are thought to be embedded in 
a large bubble of hot, low density gas (T « 10 6 Rf, n « 0.005 cm -3 ) that was probably 
produced by supernovae and stellar winds arising in the Scorpius-Centaurus OB association 
(Cox & Reynolds 1987; Frisch 1995). If this picture is correct, then the local interstellar 
medium (LISM) may contain an inhomogeneous mixture of clumps of older material swept 
up by the expanding bubble along with the more-recently processed material in the bubble 
itself. Thus, even in this local environment one may find regions of gas with somewhat 
different evolutionary histories and different compositions. The same processes have been 
at work throughout the history of the Galaxy; hence measurements of deuterium and 
other abundances along numerous lines of sight will help determine not only the chemical 
evolution of the Galaxy, but also the relative timescales for chemical evolution and the 
mixing of material. 

In this paper we present an analysis of the line of sight to the nearby white dwarf 
HZ43A. Companion papers will present similar analyses of the lines of sight to G 191-B2B 
(Lemoine et al. 2001), WD 0621-376 (Lehner et al. 2001), WD 1634-573 (Wood et al. 2001), 
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WD 2211-495 (Hebrard et al. 2001), BD+28 4211 (Sonneborn et al. 2001), and FeigellO 
(Friedman et al. 2001). A i overview will be provided by Moos et al. (2001). 

High-resolution spectra of the DA white dwarf HZ43A, covering the far ultraviolet 
(FUV) wavelength range 005 - 1 18 ?A were obtained with the Far Ultraviolet Spectroscopic 
Explorer (FUSE) for the purpose of studying the deuterium abundance of the local 
interstellar medium (ISM). This line of sight is promising for a study of D/H for several 
reasons: the star itself is moderately bright, exhibits a nearly featureless continuum, and 
has been accurately modelled; the absorbing interstellar gas appears to have only a single 
velocity component and the column density is very low, so many absorption lines that aie 
ordinarily on the flat part of the curve of growth are either unsaturated or exhibit only 
mild saturation effects. In addition, the Hi column density along the line of sight can be 
determined by two independent methods: from the shape of the Lyman continuum observed 
by the Extreme UltravioUt Explorer (EUVE), and from the Lyman-a profile observed with 
the Goddard High Resolution Spectrograph (GHRS) onboard the Hubble Space Telescope 
(HST). Systematic uncertainties in determination of Hi column densities are often the 
limiting factor in measurements of the D/H ratio, so the availability of multiple independent 
measurements of N(H) is quite valuable. 

In Section 2 we discuss the line of sight to HZ 43 A and the properties of the star; 
in Section 3 we describe the observations and data reduction procedures; in Section 4 we 
describe our analysis procedures, our methods for minimizing systematic errors, and our 
measured column densities; and in Section 5 we summarize our results. 
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2. Line of Sight and Stellar Properties 

The hot DA white dwarf HZ43A (WD 1314+293) is a member of a non-interacting 
binary system with a dMe star. The cool secondary star has no significant UV flux. The 
coordinates of the star and some other basic characteristics are given in Table 1. The 
Galactic coordinates, l = 54.10°, b = +84. 16°, are close to the north Galactic pole and the 
line of sight exits the LIC after traversing a very short distance: ~0.1pc (Redfield & Linsky 
2000). The mean H I density in the LIC is 0.1 cm -3 (Linsky et al. 2000), so the corresponding 
column density is only ns 3.1 x 10 16 cm~ 2 . The velocity of the LIC projected onto the 
HZ43A line of sight is -8.9 km s -1 (all velocities given in this paper are heliocentric). The 
line of sight may intersect the “G” cloud that lies between us and the Galactic center, and 
the “NGP” cloud that lies between us and the North Galactic Pole (see Sfeir et al. 1999 and 
the review by Frisch 1995 for a discussion of the geometry of the nearby clouds lying within 
the local bubble). The boundaries of the clouds are not known in detail, but the identity of 
each absorber can be determined by its relative velocity. In particular, the velocity of the 
G cloud projected onto the HZ43A line of sight is -12.1km s -1 . The distance to HZ43A 
has been determined by several methods to be approximately 65pc. Vennes et al. (1997) 
obtained 67pc from a comparison of its apparent and absolute V magnitude (derived from a 
model atmosphere), while Dupuis et al. (1998) obtained a somewhat smaller value, 56-61pc 
using a similar method. Ground-based parallax measurements have given 48-91pc (Margon 
et al. 1976), 53-83pc (Dahn et al. 1982), and 55-81pc (van Altena et al. 1995). The only 
discrepant measurement, 25 - 44pc, is from the Hipparcos catalog (Clark & Dolan 1999); 
the origin of this discrepancy is not understood. The apparent line-of-sight velocity of the 
star (including gravitational redshift) was measured by (Reid 1996) to be +20.6 km s -1 . 

The atmospheric parameters of HZ 43A are difficult to measure at visible wavelengths, 
because of the presence of the bright M star companion only ~ 3" away. Napiwotzki et al. 



(1993, N93 hereafter) and Finley et al. (1997, FKB hereafter) obtained optical spectra 
of both components and corrected the observed flux of the white dwarf by subtracting 
the contribution from the M star. N93 obtained T e f f = 49,000 K and logy = 7.7 using 
NLTE models, while FKB determined T eS = 50,800 K and logy = 7.99 using LTE models. 
According to Napiwotzki et al. (1999), the difference between both studies may illustrate the 
correction in T^ff and log y that one must apply to transform the atmospheric parameters 
determined by LTE models to NLTE models. However, if we apply the Napiwotzki et al. 
(1999) correction (A logy = -0.13) to FKB’s value (7.99), we obtain a larger gravity than 
N93 (7.7). 

Other spectral bands have been used to evaluate the effective temperature and gravity 
of HZ43A. For example, Holberg et al. (1986) fitted the /t/ELyman-a line profile, obtaining 
a larger temperature and gravity than the optical observations performed by N93 and FKB. 
Both Vennes (1992) and FKB indicated that the visual magnitude estimated by Holberg 
et al. (1986) was too faint. The visual magnitude is quite difficult to determine from the 
ground because of contamination by the companion star. If Holberg et al. (1986) had used 
the magnitude determined from HST data by Bohlin et al. (1995), the effective temperature 
would have been T eff = 50, 850 K, very close to that determined later by FKB. The complete 
Lyman series of HZ43A was analyzed by Dupuis et al. (1998) using ORFEUS and by Kruk 
et al. (1997) with the Aslro-1 Hopkins Ultraviolet Telescope, while the EUV wavelength 
range was studied by Dupuis et al. (1995) and Barstow et al. (1995) using EUVE. If we 
take the average of the eifective temperature and gravity for all measurements made from 
EUV, FUV, and optical observations, we obtain T e s = 50,300 K and logy = 7.88, with 
standard deviations of 700 K and 0.15 dex (see, e.g., Table 1 of Dupuis et al. 1998). 

For the purpose of this study, we adopted the following atmospheric parameters. 

T eff = 50,900 K, and logy = 8.0. These are the pre-publication values of FKB used by 
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Kruk et al. (1999) in the final Astro-2 calibration of the Hopkins Ultraviolet Telescope , and 
used presently in the FUSE calibration (the present FUSE flux calibration is defined by 
HZ43A, G 191— B2B, and GD246). 

HZ 43A is unusual in that its photosphere appears to consist entirely of hydrogen, 
though stars with similar atmospheric parameters show traces of heavy elements. Barstow 
et al. (1995) set an upper limit on the helium abundance of 10~ 3 * * * 7 , based on the EUVE 
spectrum of this star, and Dupuis et al. (1998) used an ORFEUS spectrum to set abundance 
limits of 1CT 7 for S, 10“ 8 for C and N, and 1(T 8 - 5 for Si, P, and Cl. Vennes et al. (1991), 
N93, and Dupuis et al. (1995) came to the same conclusion, though with somewhat higher 
upper limits. Therefore, photospheric models for this star employing a pure hydrogen 
composition should give accurate results. 


3. Observations 

3.1. FUSE Observations 

FUSE comprises four independent co-aligned telescopes and spectrographs. Two of the 
channels contain optics coated with aluminum and overcoated with LiF, and the other two 
channels contain optics coated with SiC. Together the four channels span the wavelength 
range 905-1187A. Most wavelengths in this band are sampled by at least 2 channels, and 
roughly one-third of the bandpass is sampled by 3 or 4 channels (1000A < A < 1080A). 

Comparison of features in the spectra obtained in different channels can thus be used to 
test for the effects of fixed-pattern noise, backgrounds, or other instrumental characteristics 

that vary from one channel to another. Each spectrograph has 3 entrance apertures: 

LWRS (30"x 30"), MDRS (4"x 20"), and HIRS (1.25"x 20"). Each aperture illuminates 
different regions of the detector and produces spectra with slightly different line-spread 
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functions (LSFs). Thus, if a source is observed through multiple apertures one can perform 
additional tests for the piesence of fixed-pattern noise introduced by the detector, and one 
can determine if contamination by geocoronal emission is significant. Further information 
on the FUSE satellite and instrument on-orbit characteristics are provided by Moos et al. 
(2000) and Sahnow et al. (2000). 

Three FUSE observations of HZ43A have been obtained, two in the LWRS aperture 
and one in the MDRS aperture. Basic information for the datasets is given in Table 
2. The data were obtain d in “histogram” mode, in which a two-dimensional spectral 
image is accumulated and downlinked. The LWRS observations will be described here for 
the sake of completeness: the signal-to-noise of the LWRS observations was too low in 
comparison with the MDRS data for them to significantly constrain the results. They were 
useful primarily for tests of fixed-pattern noise effects. Exposure durations for the LWRS 
observations ranged from just under 500 s to just over 800 s. The exposures in the MDRS 
observation were all appi oximately 485 s, apart from a few exposures that were shortened 
to roughly 200 s by delayed target acquisitions. The MDRS observation was split into 4 
intervals of approximately equal exposure time, each with a different X-position of the 
Focal Plane Assemblies (FPAs). The shift in FPA position causes a corresponding shift in 
the position of the spectium on the detector along the dispersion direction. The 4 positions 
were at X = -57 pm, 0 pm, +33 pm, and +125 pm with respect to the nominal position, 
corresponding to spectral shifts of roughly -12, 0, +5, and +19 pixels. These shifts span a 
range several times larger than the width of typical interstellar absorption lines, and several 
times larger than the width of the high-frequency components of the detector fixed-pattern 
noise. This so-called “focal-plane split” (“FP-split”) procedure results in a considerable 
improvement in the sign.d-to-noise ratio achieved, even if the the spectra from the different 
FPA positions are simply shifted and added with no attempt to explicitly solve for the 


detector fixed-pattern noise. 
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The data were processed using the most recent version of the FUSE calibration 
pipeline, CALFUSE version 1.8.7. Because there is no information available for the arrival 
time of individual photons in histogram mode, the pipeline applies exposure-averaged 
corrections for the Doppler shift induced by the spacecraft orbital motion and for the 
rotation of the gratings on orbital timescales. The primary mirrors also rotate slightly on 
orbital timescales, which can lead to significant changes in the position of the image of the 
source at the spectrograph entrance aperture. As a result, the effective spectral resolution 
can vary from one exposure to another and the zero-point of the wavelength scale varies 
from one exposure to another. 

The rotations of the primary mirrors during an exposure can cause the stellar image 
to fall outside of the corresponding spectrograph aperture, resulting in loss of signal in that 
channel for all or part of an exposure. The Fine Error Sensor presently used to provide 
pointing information to the satellite attitude control system is part of the optical path in 
the LiFl channel; hence the target is always maintained at a fixed position in the LiFl 
spectrograph aperture. The alignment of the other channels is usually maintained so that 
the rotations of the primary mirrors occurring on orbital timescales do not shift the source 
image at each focal plane out of the corresponding LWRS aperture. For the M1010501 
observation, no flux was lost in any of exposures; for P1042301 no flux was lost except 
that the SiCl channel fluxes were low by 7%, 10%, and 19% in exposures 13, 16, and 
20, respectively. The SiCl spectra obtained in these exposures were scaled accordingly 
before being combined with the other exposures. For MDRS observations, peakups are 
performed twice in each orbit in order to try to correct for the mirror rotations by moving 
the spectrograph apertures (FPA)s in the dispersion direction to maximize the signal prior 
to the start of an exposure. For the P1042302 MDRS observation, this procedure was fairly 
successful, providing effective exposure times of 82% for LiF2, 62% for SiCl, and 76% for 
SiC2. The necessity of moving the FPAs to compensate for rotations of the primary mirrors 
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means that there is actually a distribution of FPA positions clustered about the mean 
offsets that define the FP split sequence described above. 

The spectra obtained in different exposures generally need to be co-aligned (shifted 
to a common wavelength zero-point) prior to combining. This is caused in part by the 
lack of corrections for primary mirror rotations in the present pipeline, and in part by 
shortcomings in the prese nt corrections for grating rotations. For HZ 43 A there are no 
narrow photospheric lines and very few strong interstellar absorption lines. For the LiF 
channels, we fit a Gaussian profile to the interstellar C II 1036. 33?A line in each exposure 
and shifted each exposure to locate this line at a common velocity. For the SiC channels, we 
performed a cross-correlation on the the interstellar Lyman lines between 915A and 930A 
(fitting the Cll 1036.337A line gave essentially identical results). In order to have enough 
signal to measure line positions, we discarded spectra for which the flux was below half of 
the nominal value for its channel. 


3.2. GHRS Observations 

Archival tfST/GHRS spectra were used to determine the velocity structure of the 
intervening gas and to measure Lyman-c* absorption of both Hi and Dl. A log of the GHRS 
observations is given in liable 3. All of the data were acquired after the installation of the 
Corrective Optics Space Telescope Axial Replacement (COSTAR) unit. The echelle-mode 
data, all of which were acquired through the large science aperture (LSA; 1.74" x 1.74"), 
have a line spread function with a narrow core of resolution 3.5 km s -1 (FWHM) and broad 
wings (-15 km s _1 FWHM) with relatively low power (see Appendix B of Howk et al. 

1999). 


Our reduction of these data follows that of Howk et al. (1999). The basic calibration 
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makes use of the standard CALHRS pipeline with the best calibration reference files 
available as of the end of the GHRS mission. The CALHRS processing includes conversion 
of raw counts to count rates and corrections for particle radiation contamination, dark 
counts, known diode nonuniformities, paired pulse events and scattered light. The 
wavelength calibration was derived from the standard calibration tables based upon the 
recorded grating carrousel position. The absolute wavelength scale should be accurate to 
~1 resolution element, or ~3.5 km s -1 for the echelle-mode data used in this work. 

All of the echelle-mode observations employed here were taken at four separate grating 
carrousel positions to place the spectrum at different positions along the 500-diode science 
array of the GHRS. This is done to mitigate the effects of the fixed-pattern noise that can 
limit the signal-to-noise ratio and produce artifacts in GHRS spectra. We have co-added 
each of the individual sub-exposures by aligning the interstellar lines of interest (e.g., Dl, 
Mg II, and Fell) in each of the observations. We used the algorithm of Cardelli & Ebbets 
(1994) to explicitly solve for the fixed-pattern noise spectrum in each of the datasets. 
However, we found it unnecessary to apply this correction because the signal-to-noise ratio 
of the data did not necessitate the removal of the low-power features and no prominent 
features were found (see discussion in Howk et al. 2000). We note that the coaddition of 
the shifted datasets reduces the strength of any fixed- pattern features by a factor of four. 

The inter-order scattered light removal in GHRS echelle-mode data discussed by 
Cardelli et al. (1990, 1993) is based upon extensive pre-flight and in-orbit analysis of GHRS 
data and is used by the CALHRS routine. The coefficients derived by these authors are 
appropriate for observations made through the small science aperture (SSA). Corrections to 
these coefficients are often required to bring the cores of strongly saturated lines observed 
through the LSA to the correct zero level. Typically this is done by adjusting the spectrum 
zero level by an amount equal to a fraction of the average net flux in each order. This “d” 
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coefficient (Cardelli et al. 1990, 1993) accounts for the residual effects of the local echelle 
scattering. 

For the Ech-B data covering Mgn and Fell, no strongly saturated lines were available 
for evaluating the quality of the background removal, so we adopted the standard 
d-coefficients from the CAL HRS processing. The Lyman-a absorption profile for HZ43A 
should be flat-bottomed at zero flux near line center. Although the presence of strong 
Lyman-a airglow confuse-; the matter somewhat, it is clear that the use of the standard 
d-coefficient (d=0.052) over-subtracts the core of the Lyman-a profile. We have used 
the observed flux level over the wavelength range 1215.74 to 1215.83A, assuming this 
represents the true zero point, to adjust the d-coefficient. The adopted d-coefficient is d = 
-0.006i0.033. 

The G160M observations of the 1192-1228A range were acquired through the Small 
Science Aperture (0.22"x0.22"; SSA). These data have a resolution of -19.6km s' 1 
(FWHM). The scattered light properties of this holographically-ruled first-order grating 
are excellent. No corrections to the CALHRS background subtractions were necessary. 


4. Analysis 
4.1. Method 

We attempted to determine column densities by more than one method whenever 
possible. The two primary methods were profile fitting and curve-of-growth analyses. 
These two methods differ in their sensitivity to some of the sytematic uncertainties that 
affect column density measurements. A curve-of-growth analysis was not feasible for Hi, as 
essentially all of the lines other than Lyman-a are on the flat part of the curve of growth. 
Instead, H I column density measurements based on fits to the EUV spectrum of HZ 43A 
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were taken from the literature to provide a consistency check on our results. 

Profile fitting was used to determine the column densities of Hi, Dl, O I, and Nl. The 
profile fitting was performed using the program Owens. f, written by one of us (M.L.). This 
program uses x 2 minimization to adjust temperatures, turbulent velocities, and relative 
velocities for multiple gas components, to adjust column densities for each species present 
in each gas component, and to adjust polynomial fits to continuum fluxes, line spread 
functions, background levels, etc. for each specified spectral region. Each quantity can be 
adjusted by the program or held fixed to user-specified values. Further information on this 
program is provided in the companion paper by Hebrard et al. (2001). 

The FUSE and GHRS spectra were normalized prior to fitting by means of an NLTE 
stellar model. The continuum in the vicinity of individual absorption lines was then 
fit by a low-order polynomial as required to account for any residual flux calibration 
errors or flat-field effects. We use the program TLUSTY (Hubeny &: Lanz 1995) to 
calculate the atmospheric structure of HZ43A under the assumption of NLTE, and then we 
compute the detailed emergent flux using the program SYNSPEC (Hubeny 2000, private 
communication). The hydrogen line profiles are obtained from the Stark broadening tables 
computed by Lemke (1997) using the program of Vidal et al. (1973). We use a NLTE model 
atmosphere to calculate the emergent flux, because NLTE effects are important in the core 
of the Lya line, which is formed high in the atmosphere where departure from LTE is 
significant (see, e.g., Wesemael et al. 1980). The model parameters used were as discussed 
in Section 2: the composition was pure hydrogen, T e ff= 50900K, and logg= 8.00. The 
model spectrum was then normalized to V = 12.914 (Bohlin et al. 1995). Except for NLTE 
effects in the narrow cores of the Lyman lines it is essentially identical to the LTE model 
described in detail by Kruk et al. (1999) and references therein. 

The Lyman-a profile was also fit using a different program written by one of us 
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(J.C.H), and the profile fitting code of Ed Fitzpatrick (Spitzer k Fitzpatrick 1993) was used 
to provide a consistency check for certain test cases. 

We have derived column densities for N I and O I by fitting the measured equivalent 
widths with a single-component Gaussian (Maxwellian) curve of growth (Spitzer 1978, 
Jenkins 1986). The analysis was performed using software adapted from that described 
by Savage et al. (1990), which varies the Doppler parameter, b, and column density so as 
to minimize x 2 between the observed equivalent widths and a model curve of growth. The 
stellar continuum in the vicinity of each line was estimated using a low-order Legendre 
polynomial fit to the data. The dispersion of the data about the fit was used to estimate 
the uncertainties in the final integrated equivalent width measurements as discussed by 
Sembach k Savage (1992). This estimate includes contributions from both Poisson noise and 
uncorrected high-frequen y fixed-pattern detector features. The measurement uncertainties 
quoted include contribut ons from this noise estimate, uncertainties in the Legendre fit 
parameters (Sembach k Savage 1992), and estimates of the systematic uncertainties 
inherent in the continuum placement and velocity integration range. Atomic data for both 
the profile fitting and the curve of growth analysis were taken from Morton (1991, Morton 
2001, private communication). 

The integrated equivalent widths for several interstellar species are presented in Table 
4, and the results from the curve-of-growth fits are given in Table 5. For most species other 
than N I and O I the number of available transitions was insufficient to fit a curve of growth 
and determine a b- value; instead, the measured equivalent width was combined with a curve 
of growth computed using the O I b- value to derive the column density. For these cases 
the table also lists the difference between the column density derived in this manner and 
that derived using a straight integration of the apparent column density profiles (Savage k 
Sembach 1991), providing a measure of the saturation correction. 
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4.2. Systematic Uncertainties 

Column density determinations are subject to numerous possible systematic 
uncertainties. These are described in the following section, along with brief descriptions of 
how these were addressed. The magnitude of each effect on the column densities will then 
be addressed in the discussion of each species as appropriate. 


4-2.1. Overlapping velocity components along the line of sight 

One of the major difficulties typically encountered in measuring column densities is 
accounting for overlapping absorption profiles resulting from clouds with similar velocities. 
As we will show below, there appears to be only a single absorption component along the 
line of sight to HZ 43A, so this complicating factor is not present. 


4-2.2. Narrow absorption features in the stellar spectrum 

Stellar photospheres ordinarily contain trace quantities of numerous elements that can 
give rise to narrow absorption features in the emergent spectrum. In some cases the stellar 
spectrum is well-understood and the presence of such narrow features can be predicted 
and possibly even modelled. In other cases, the wavelengths of the various transitions 
are poorly known and it is impossible to be certain how much of an absorption profile is 
due to interstellar gas and how much results from the photospheric spectrum. However, 
as discussed in Section 2, the atmosphere of HZ43A appears to be composed solely of 
hydrogen. Thus, this potential complicating factor is also not a concern for our analysis of 
HZ43A. 
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4-2.3. Uncertainties in the stellar Lyman-oc profile 

The high gravity of HZ43A causes the photospheric Lyman lines to be quite broad, 
hence the stellar spectrum is varying slowly in the vicinity of all of the interstellar lines 
being fit other than H I Lyman-a. Because we are fitting the damping wings of the 
interstellar H I Lyman-a profile, spanning roughly 6A, we need to be able to model the 
photospheric absorption very accurately. In addition, the H I column density along the line 
of sight to HZ43A is sufficiently low, and the relative velocity of the star and interstellar 
gas is sufficiently high, that the NLTE core of the photospheric Lyman-a profile may affect 
the red wing of the combined profile. HZ43A has been studied in great detail, however, and 
model atmospheres and synthetic spectra of HZ43A, along with those of several similarly 
well-studied DA white dwarfs, are used to define the flux calibration for most orbiting UV 
spectrometers (Bohlin 2000, and references therein). The parameters describing HZ43A are 
discussed in Section 2, and the fits to the Lyman-a profile are discussed in Section 4.3.2. 


4-2-4- Continuum placement errors 

Errors in placement of the continuum or the background will clearly affect the column 
density determination. In most regions of the spectrum the continuum appears to be smooth 
and slowly-varying in the vicinity of each of the lines measured. If the continuum did not 
appear to be well-behaved, most likely because of large-scale fixed-pattern noise effects 
in the detector, then the the affected line was not included in the fit (the SiClB MDRS 
measurement of the O I D48.686A line is such a case). Uncertainties in continuum placement 
are included explicitly in the equivalent width measurements used in the COG analysis, and 
implicitly in the x 2 minimization performed by the profile fitting. This latter approach is 
adequate if the true continuum is in fact well-represented by a low-order polynomial, and 
if the absorption being measured occupies a small fraction of the region being fit. For the 
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broad Lyman-a profile a variety of continuum fits were tested to determine the sensitivity 
of the results to the adopted continuum. This will be discussed in greater detail below r in 
the section on the measurement of the H I column density. 


4-2.5. Background placement errors 

The dominant background for FUSE spectra of a bright source such as HZ 43A is 
grating scattering of source flux. The intrinsic detector background and diffuse stray light 
in the spectrograph are negligible in comparison. The grating-scattered source flux has 
been characterized by measuring the flux in the cores of highly-saturated H 2 absorption 
line profiles in the spectra of stars with translucent clouds along the line of sight. Typical 
results for this background are that the level is about 1% of the continuum flux in the SiC 
channels, and 0.3% to 0.5% of the continuum flux in the LiF channels. In the HZ43A 
FUSE observations, the flux shortward of the Lyman limit is indeed about 1 % of the 
continuum flux in the SiClB channel. None of the interstellar absorption lines in the 
HZ43A spectrum are broad enough for the true background level to be seen in the line 
cores. The instrumental LSF has a broad component that has a FWHM of about 18 pixels, 
which is broad enough that even the H I Lyman line cores do not reach the level of the 
background. This point is discussed further in Section 4.2.8. The effects of uncertainty in 
the background level were investigated by performing the profile fits with the background 
level set to zero, set to our best estimate of the background, and and set to twice our best 
estimate of the background. The resulting column density variation was typically 0.01 
dex, and is included in the quoted uncertainties for the column densities. The background 
effects in the GHRS spectra are discussed in the section below on the Hi column density 


measurement. 
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4-2.6. Geocoronal emission 

Contamination of the spectrum by geocoronal emission can be significant for the large 
spectrograph apertures employed by FUSE. Such emission is quite strong at Lyman-/?, but 
declines rapidly in strength for the higher Lyman lines. Geocoronal emission is also present 
on the day side of the orbit for many of the same O I and N I transitions for which we are 
trying to measure interstellar absorption. The widths of geocoronal lines are limited by the 
size of the spectrograph apertures. For the LWRS aperture this width is typically 0.35A, or 
the equivalent of 106 km s' 1 at Lyman-/?, while for the MDRS aperture this width is about 
0.045 A, or the equivalent of 13 km s _1 at Lyman-/?. 

The MDRS aperture is sufficiently narrow (~13km s -1 ) that the geocentric velocity of 
the interstellar gas shifts the geocoronal emission almost completely away from narrow lines. 
The red wing of the broad interstellar Lyman-/? line is filled in by geocoronal emission, but 
no difference is seen in the Ol or Ni absorption profiles when the day and night data are 
analyzed separately. The absorption profiles measured in the LiFl channel are actually 
deeper for the day data than the night data, presumably because the LiFl grating is more 
stable in the day than in the night and the effective resolution is correspondingly better. 


4.2.7. Uncorrected fixed-pattern noise 

Uncorrected fixed-pattern noise on small scales can lead to significant errors in the 
shapes and fluxes of line profiles; on larger scales it can lead to continuum placement 
errors. The FUSE detectors exhibit significant fixed-pattern noise on both small and large 
scales, caused primarily by the inherent granularity of the microchannel plates (MCPs). 
The slowly-varying relative spacing of the pores in the stack of 3 MCPs gives rise to Moire 
fringe effects with a spacing of 8-9 pixels, comparable to the width of the LSF. Changes in 
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pore spacing at fiber bundles and regions of dead pores can cause changes in the detector 
quantum efficiency on larger scales. The present CALFUSE pipeline flags pixels affected by 
dead regions of the detector, but otherwise fixed-pattern features are not yet corrected. 

We have mitigated the effects of fixed-pattern noise by several means. Every transition 
being examined in the FUSE bandpass is sampled by either 2 or 4 channels, so we have 
verified that each such measurement provided essentially similar results. In addition, 
because HZ 43A was observed in both LWRS and MDRS apertures, we have checked for 
consistency of the results. The MDRS data were obtained with an FP-split procedure, 
which substantially reduces the fixed-pattern noise on scales of one to two resolution 
elements. The resulting spectra typically have a S /N that is about 80-90% of that expected 
from photon statistics, depending on the size of the region examined, primarily because the 
larger-scale fixed-pattern noise has not been removed. However, most of these fluctuations 
are slowly varying and can be removed by low-order polynomial fits to the continuum in 
the vicinity of absorption lines of interest. 


4-2.8. Uncertainties in the instrumental line spread function 

Characterization of the FUSE LSF is still in a preliminary state. There appears to be 
both a narrow and a broad component to the LSF. The narrow component is somewhat 
wavelength-dependent and is sensitive to the accuracy of the correction for grating motion 
during exposures and to the accuracy with which the user was able to co-align the spectra 
from each exposure prior to combining them. In addition, data acquired in histogram mode 
cannot be corrected for the degradation of the spectral resolution caused by rotation of 
the primary mirrors or gratings during an exposure. For example, in some channels the 
resolution is better during the night portion of the orbit and for other channels it may be 
better during the day. As a result, the effective width of the narrow component of the LSF 
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will vary from one observation to another. We measured the effective width of this narrow 
component by fixing the temperature and turbulent velocity of the interstellar gas to that 
found from the GHRS measurements (see below), and letting the LSF widths vary in the 
profile fits. The results were generally consistent from one absorption line to another, and 
were typically about 9 pixels FWHM (corresponding typically to 0.062 A in the LiF channels 
and 0.057A in the SiC channels). This is consistent with what is found from other FUSE 
datasets. These widths were then held fixed when performing subsequent fits. It should be 
noted that the LSF width is significantly greater near the edges of the detector active area, 
as are the uncertainties in the wavelength calibration. 

The broad component of the LSF appears to be somewhat less sensitive to the 
effects of grating rotation, primary mirror rotation, or coalignment of spectra. The initial 
characterization the available data indicates that about 30 percent of the area of the LSF 
is in a component with a width in the range of 17-24 pixels. Analysis of highly saturated 
lines that are flat-bottomed over regions of 150 pixels or more reveals that there is no 
significant component of the LSF with a characteristic width larger than 25 pixels. The 
effects of this broad component of the LSF on the profiles of unsaturated lines may not be 
evident unless the signal- to-noise ratio of the data is very high; however, its main effect is to 
introduce a substantial apparent increase in the background under the line. If a weak line 
is fit with a single, narrow, Gaussian LSF, the background level in the fit must be increased 
to account for this effect In the profile-fitting analysis of the HZ 43A ISM absorption we 
have represented the LSI as a two-component Gaussian function. The width of the narrow 
component was determined to be about 9 pixels as described in the previous paragraph, 
and the broad component width was set to 17 pixels, with an amplitude corresponding to 
30% of the total area. 

If the LSF has very long shallow tails that cannot be distinguished from the continuum, 
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then both the profile fitting technique and curve-of-growth analysis would likely be subject 
to continuum placement errors that will cause some of the area of the line to be missed and 
therefore cause the the column density to be underestimated. However, while there is a 
component of the LSF that has roughly twice the width of the narrow component, there is 
no component with very long extended tails. In the absence of such a weak and very broad 
component, the curve-of-growth analysis is independent of the shape of the LSF, and thus 
provides a valuable check that the profile fit was not adversely affected by errors in the 
assumed LSF. 

The LSF of the GHRS modes used here is smaller than the Doppler widths of the 
lines being measured, other than for Fell, so there is no significant systematic error 
associated with uncertainties in this LSF. The Doppler width of the Fell 2382. 765A line is 
comparable to the LSF width, so the determination of the Doppler width for this line does 
have relatively large uncertainties. The column density, however, is not adversely affected 
because the line is so weak. 


4-2.9. Uncertainties in the f-values 

Errors in the atomic data used in the fits have the potential of biasing the results. 

If such errors are distributed randomly then the use of numerous lines will minimize any 
systematic bias of the overall column density determination. If the oscillator strengths for 
a given species are all too large or too small then the estimate for the column density will 
be biased. However, such a systematic error in the f-values seems unlikely to result in a 
consistent fit to all of the lines, especially if lines are included that are not on the linear part 
of the curve of growth. We have tested for effects on the derived column density resulting 
from errors in the f-values by repeating the profile fits after removing each transition and 
verifying that the derived column density did not vary by more than 1 a. An examination of 



- 23 - 


the curve-of-growth plots also shows that a satisfactory fit can be obtained that is consistent 
with each of the included points. The only transition that was found to have an f-value 
clearly inconsistent with those of the other transitions was the that of the Ol 1026. 473A 
line (see Sembach 2000). This line was excluded from all fits and from the curve-of-growth 
analysis. 


4.3. Results 

4.3.1. Line of sight velocity structure 

The velocity structure of absorbing gas along the line of sight was determined from 
GHRS observations of Dl 1215. 339A, Fell 2383.765A, Mg II 2796.352A, and Mgll 2803. 53lA 
line profiles. These lines are w r ell-resolved, and are isolated from any other absorption 
features. We obtained a satisfactory fit to the lines with a single absorption component 
with a heliocentric velocity of -5.2 ± 0.1km s -1 , a temperature of T = 5353 ± 948 K, 
and a non-thermal velocity of h nt = 2.0 ± 0.2km s -1 (the uncertainties are lcr statistical 
uncertainties; the systematic uncertainty on the absolute velocity is ± 3 km s x )- The total 
Doppler parameter, b, is related to the temperature and non-thermal velocity parameter 
such that b 2 = 0.0165 Tf l + b 2 nt . The observed profiles and the best model fits are shown 
in Figure 1. There is no evidence for an additional velocity component. These values for 
T and b nt were not held fixed in subsequent fits (e.g. for Ol or Nl). Instead, the derived 
Doppler widths for each species were computed and compared with the ranges allowed by 
the above values; in each case the derived Doppler widths were found to be consistent with 
the la uncertainties quoted above. The b- values for Ol and Nl predicted from the values 
of T and b nt given above are: b(Ol) = 3.1 ± 0.2km s *and b(Nl) = 3.2 ± 0.2km s {la 
uncertainties) . 



-24- 


4-3.2. Hi Column Density 

The interstellar H I column density along the line of sight to HZ 43A can be determined 
by three different means: fitting the spectral shape of the EUV continuum, fitting the 
damping wings of the H I Lyman-a profile, and fitting the converging Lyman series. This 
last method has the greatest uncertainties because the lines are typically all saturated, but 
for the low column density of this line of sight it is possible to obtain meaningful results. 

Determination of the H I column density by analysis of the EUV spectrum of a star 
requires accurate knowledge of the stellar continuum flux. The EUV spectrum of a star is 
very sensitive to the presence of even trace quantitites of heavy elements in the photosphere, 
and substantial effort is usually required in order to model this spectrum accurately. HZ 43A 
is unusual in that no traces of heavy elements have yet been detected in its photosphere; 
hence it is not affected by this potential source of systematic error. EUVE spectra of 
HZ43A have been analyzed by several groups, using different observations, different model 
atmosphere codes, and different fitting methods. The three most recent analyses of EUVE 
observations will be described briefly below. 

Dupuis et al. (1995) analyzed observations from 1994 March 25,28 and 1994 May 20. 
They first determined the He I and He II column densities by fitting the ionization edges 
at 504A and 228A, respectively, and then determined the stellar T e g- and interstellar H I 
column density by fits over the 400-550A bandpass. This wavelength range was chosen 
because it was least sensitive to variations in surface gravity or heavy element abundance. 
The fitting was done with model spectra computed from pure H model atmospheres 
(Vennes & Fontaine 1992). The results were T e ff = 51100K ± 500K, and Nhi = 

(8.7 ± 0.6) x 10 17 cm~ 2 (1 a). They also estimated that systematic uncertainties in the 
effective area of 20% would change T e g by 3% and Nhi by 10%. Uncertainties in the 
V magnitude (used to normalize the models) would have similar effects, but are small in 
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comparison. The ionization fraction of hydrogen was estimated by assuming that the ratio 
of total H to total He was 10:1, that there was no appreciable He III, and comparing the 
measured N H i to the measured He I and Hell column densities. The resulting ionization 
fraction for H, fn = N(H il) / N(H), was < 0.4 ±0.1. 

Barstow et al. (1997 analyzed the 1994 May 20 observation also, but used H+He 
models computed by Koester (1991) and determined all stellar and ISM parameters from 
a simultaneous fit to the entire EUVE spectrum. The values for T e ff and logg were 
constrained to remain within the 1 a uncertainties of the values determined by Napiwotzki 
et al. (1993). The V magnitude used, 12.99, corresponds to model fluxes that are fainter 
by a factor of 1.077 than would be the case if the present best value of V = 12.909 were 
used (Bohlin 2000). This would result in an underestimate of the derived N H i column 
density by a few percent. They performed fits with two classes of models, homogeneous 
H+He mixtures, and stratified compositions. The derived interstellar Hi column densities 
and 1 a uncertainties were Nhi = (8.8 ± 0.2) x 10 1 'cm 2 (homogeneous) and Nhi — 
8.3loi x 10 17 cm -2 (stratified). The difference between the results for the two models 
exceeds the statistical errors in the fits, indicating that systematic uncertainties in the 
nature of the appropriate model are significant. The ionization fraction w r as derived in the 
same manner as by Dupuis et al. (1995) to be fn — 0.19^ Jg for homogeneous models, 
and f H = 0. 22^0^4 for stratified models. We have assumed that the results presented by 
Barstow et al. (1997) supercede the earlier results in Barstow et al. (1995), so we will not 
discuss the earlier results here. 

Wolff et al. (1999) analyzed observations obtained 1997 June 25,26. They first fit 
the He I and Hell ionizal ion edges to determine the corresponding column densities, and 
then fit the full spectrum to determine the H I column density. The stellar spectrum was 
computed using the LTE model atmosphere code of Koester (1996), assuming a pure 
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hydrogen composition. The atmosphere parameters were not free, but fixed at values 
derived from fits to the Balmer lines. The values used for T e g were 50800 ± 300 K, and the 
corresponding column densities were N H i = 8 . 851 ^ 03 ^ xl0 17 cm -2 . The uncertainties in 
the column densities are dominated by the uncertainty in the effective temperature (Wolff 
2001, private communication). The V magnitude used to normalize the models was 12.89; 
this would overestimate the stellar flux by about 2% and would result in only a very slight 
overestimate of the H I column density. The ionization fraction was derived in the same 
manner as the previous two analyses to be fn = 0.031q!o3- 

An EUV spectrum of this star was also obtained by HUT and analyzed by Kimble 
et al. (1993) to derive N H i = (6.5 ± 0.5) x 10 17 cm~ 2 (lcr statistical errors; systematic 
uncertainties were estimated at 15-20% and are dominated by calibration uncertainties). 
This falls below the more recent EUVE determinations by approximately the combined 
random and systematic uncertainties. These uncertainties are considerably larger than those 
of the much higher signal-to-noise EUVE spectra, so we will not consider this measurement 
further. The three EUVE measurements described above will be addressed again at the 
end of this section. 

We fit the GHRS spectrum of the interstellar H I plus D I Lyman-a profile towards 
HZ 43A to determine the properties of the neutral hydrogen in this direction. The spectrum 
was first scaled by a factor of 0.866 in order to match our model flux over the regions 
1212.25 - 1213 . 25 A and 1218.0 - 1218 . 5A. The spectrum was then normalized by an NLTE 
stellar model prior to fitting (see Section 4.1 for a description of the model). We found that 
additional normalization by a low-order polynomial was required to account for potential 
small wavelength-dependent errors in the flux calibration. The apparent photospheric radial 
velocity of HZ43A ( v ph ) was measured by Reid (1996) to be +21 km s -1 . The spectrum and 
the model (shifted by +21 km s” 1 ) are shown plotted in Figure 2. A visual inspection of the 
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figure shows that at least a linear polynomial will be required as an additional correction to 
the flux calibration. 

No explicit uncertainty was given by Reid (1996) for the radial velocity, but it is 
probably about 10-15 km s~ T We therefore examined the effects on the derived Hi column 
density of varying this velocity. We created stellar models with the photospheric velocity 
set at zero, at the nominal value of +21km s -1 , and at ±5, ±10, ±15 km s 1 with respect 
to the nominal value. The GHRS Lyman-d spectrum was normalized by each model in 
turn, and the best fit to the spectrum was generated for fixed values of log N H i ranging 
from 17.80 to 18.00 in steps of 0.01 dex. The resulting x 2 contours are plotted in Figure 3. 
The Fell and Mgll lines were fit simultaneously with the Lyman-c* Hi and D I line profiles. 
Because the Lyman-ct profile extends over several Angstroms, it may be affected more by 
uncertainties in the flux calibration than are narrow lines. To investigate this possibility, 
we generated these x 2 contours for continuum fits performed with both fourth-order and 
seventh-order polynomials. There was little correlation in either case between log Nhj and 
the photospheric velocity For the fourth-order fits, the minimum x 2 was obtained for a 
photospheric velocity of 26 km s _1 , but the x 2 increased slowly for smaller velocities and 
increased rather steeply for higher velocities. For the seventh-order continuum fits, the 
minimum x 2 was also at 26 km s — \ but it increased rather rapidly for both lower and higher 
velocities. The Hi column density and 2 o uncertainty derived was log N H i = 17.8951S 
from the fourth-order continuum fits, and log Nhj — 17.905lo.Q65 from the seventh-order 
fits. The corresponding x 2 values were 2750 for 2871 degrees of freedom, and 2711 for 2868 
degrees of freedom, respectively. Based on these fits, we have adopted the value log A hi 
17. 90122! f° r our combined GHRS and FUSE datasets. 

We examined the behavior of the continuum fits by varying the order of the polynomial 
from 1 to 8, for photospueric velocities (v ph ) of 0, ±5, ±11, ±21, and ±31 km s \ and 
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allowing the Hi column density to vary freely in each fit. For velocities of 0 and +5 km 
s -1 there was relatively little effect of polynomial order on the reduced x 2 (f° r orders > 

1). For higher velocities the x 2 improved significantly for orders up to 4 and improved only 
slightly thereafter (with minima at order 7), except that for +31 km s _1 the x 2 dropped 
again for order 8. An examination of the polynomials in each case show that the dominant 
term gradually changes from linear at low v P h to cubic at high v p h, and that the amplitudes 
increase from ± 2% at low v ph to ± 3% at higher v ph . The higher-order terms gave slight 
improvements to the overall x 2 » but did not change the character of the continuum fits. 

We evaluated the effects of the uncertainties of the background flux level on the derived 
log Nh i ■ Our best estimate of the background level at Lyman-a was 0.5% of the continuum 
flux level, obtained by letting the background be a free parameter in the fit. We then varied 
the background level by +2% and -2% of the continuum flux level, causing log N H i to 
change by only ± 0.012. 

We were unable to obtain a satisfactory fit to the observed Lyman-a profile with a 
single component of neutral material, but found that a good fit was possible only if a 
second, low-column density absorber were included. The temperature, column density, and 
relative velocity of this component were free parameters in the fits to the Lyman-a profile 
used to construct the x 2 contours described above. For the fits in which the continuum 
was represented by fourth-order polynomials, the column density of this second component 
was log Nhi = 14.9^o 5 (2 a). For the fits in which the continuum was represented by 
seventh-order polynomials, the column density of this second component was log Nhi = 
14. 9^2 (2 a); We will adopt the value log N H i = 14.9t°5 (2 a). The temperature of this 
component was 30,000 K ± 10,000 K (2 a) for both cases. As the column density of the main 
Hi component was increased, the derived column density for this second component would 
also increase, and its temperature would decrease. Similarly, as v p h was increased, shifting 
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the NLTE core of the photospheric Lyman line closer to the red side of the saturated portion 
of the interstellar absorption profile, the column density of this second component would 
increase and the temperature would decrease. The heliocentric velocity of this component 
was found to be about -2.0 km s’ 1 , but as v ph was varied from it would shift in the opposite 
direction (i.e. it would in< rease as v p h decreased). We suggest that this low-column density, 
high-temperature Hi component may be due to the “hydrogen wall about the Solar 
System. The hydrogen wall is an accumulation of gas at the boundary of the heliosphere, 
where the solar wind and LISM interact. Hydrogen walls have been detected around the Sun 
along several lines of sight, and around several other stars (Wood & Linsky 1998; Lmsky & 
Wood 1996; Wood et al. 2000). The gas in these walls is fairly hot, 30,000K - 180,000K, and 
Hi column densities range from log Nhi = 13.75 to log Nhi = 14.75. The properties of the 
second Hi component in our fits are consistent with these other detections. One potential 
problem with this interpretation, however, is that the Lyman-a profile observed for the 
similar high-latitude sightline to 31 Com (l - 115 °,6 = +89°) can be fit very well without 
any high-temperature component (Wood et al. 2000). The models of the interaction of the 
heliosphere and the LIC described by Wood et al. (2000) do predict absorption by such a 
hot component for the line of sight to 31 Com, so one possible explanation is that the LIC 
is inhomogeneous on small scales and that the corresponding conditions at the boundary of 
the heliosphere vary accordingly. 

The best fit to the normalized Lyman-a profile, including the main component and 
the high-temperature component, is shown in Figure 4. The relative velocity of the main 
absorbing component matched that found for the Dl, Mgll, and Fell lines, about -5km 
s _1 . This differs significantly from that expected for the local interstellar cloud (LIC), 
-8.8km s” 1 , and the column density is much higher than would be expected from the LIC 
along this line of sight (log Nhi — 16.47; Redfield &: Linsky (2000)). The velocity differs 
even more from that expected for the ‘G’ cloud, -12.1km s l . The velocity is consistent, 
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however, with that found for the single absorber along the line of sight to 31 Com by- 
Piskunov et al. (1997) (-4 ± 1km s _1 ). This absorber is called the North Galactic Pole 
cloud by Linsky et al. (2000, see their Figure 1). 

The FUSE spectrum of the converging Lyman series might in principle be used to 
try to determine the Hi column. Such fits give results consistent with the column density 
derived from the Lyman-a profile, but with rather large uncertainties. The main difficulty 
is that the flux calibration is inherently uncertain at precisely the wavelengths needed to 
perform the fits to the Lyman edge, because there is no reference star available for which 
unattenuated continuum flux is present at the Lyman edge. For any given set of parameters 
we find that the flux is overestimated at some wavelengths and underestimated at others, 
by about ±5%. If we allow for such uncertainties we find that variations of the gas Doppler 
b parameter within 1 a limits result in variations in the derived log Nhj of about 0.15 dex, 
which is too large to provide a significant constraint on our results. Ultimately, we may be 
able to determine a more accurate calibration at these wavelengths by means of combining 
data from several white dwarfs for which reliable H I columns have been determined by 
other means. The interstellar H I absorption at the Lyman edge predicted from the column 
densities, temperatures, and turbulent velocity derived above from the GHRS data is shown 
in Figure 5. The agreement is excellent, indicating in particular that the Doppler width of 
the gas is accurately known. The effects of the broad component of the instrument LSF is 
clear in this Figure, as the flux in the cores of the saturated Lyman lines is several times 
higher than the background level seen below the Lyman edge. 

We will combine the available EUVE - derived measurements of the Hi column density 
with our GHRS / FUSE measurement to provide our best estimate of Nhi- The EUVE 
measurement uncertainties are driven in large part by the uncertainty in the effective 
temperature. Each group used different estimates for this uncertainty, but they show 
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similar sensitivity in the sense that a change in T e ff of 500K would correspond roughly 
to a change in N H i of 0.(i x 10 17 cnT 2 (this scaling is not apparent in the Barstow et al. 
1997 paper, but in the earlier Barstow et al. 1995 paper N H i appears to be roughly 
50% more sensitive to changes in T e ff). If we take the uncertainty for T e ff from Dupuis 
et al. (1995) as representative, the corresponding uncertainty in N H i from an individual 
EUVE measurement is roughly equal to the 1 a uncertainty in our determination of N H i 
from the GHRS and FUSE data. Therefore, we will simply take the mean of the four 
measurements (after first averaging the two values obtained by Barstow et al. (199 / ) for 
both the homogeneous and stratified models), which gives N H i = 8.511 x 10 17 cm 2 , or 
log Njji = 17.930. The standard deviation of the four values was 0.398 x 10 17 , or 0.02 
dex, but this may not be an adequate representation of the true uncertainty. Each of the 
EUV measurements is affected in a similar fashion by the uncertainties in the effective 
temperature. Barstow et al. (1997) were able to obtain relatively small uncertainties for 
N H i in the context of either the homogeneous or stratified models, but the differences 
between the two corresponding column densities was 0.5 x 10 17 cm , or almost as much 

as for the nominal uncert ainty of 500 K for T e ff. It is worth noting that the two recent 
optical determinations of T eff have uncertainties greater than 500 K: 49000 ± 2000 K (3 a) 
by Napiwotzki et al. (1993) and 50822 ± 639 K (1 <r) by Finley et al. (1997). Similarly, as 
we showed above in Section 2 the standard deviation for all the recent determinations of 
Teff is about 700 K. Therefore, we prefer to use a somewhat more conservative estimate for 
the uncertainty in Nhi, and adopt the value log Nhi — 17.930 ± 0.060 (2 cr). 


4.3.3. Ol Column Density 

The FUSE bandpass contains numerous O I absorption lines of varying strength. The 
total column density to HZ43A is sufficiently low that even the strongest O I lines are close 
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to being on the linear part of the curve of growth. The data were analyzed with profile fits 
and a curve-of-growth analysis. Good quality data were obtained in all channels for the 
MDRS observation. Separate analyses of the day and night MDRS data actually give a 
somewhat higher value for N(Ol) than the night data; the opposite would be expected if 
airglow contamination were significant. 

The MDRS data were analyzed in 4 ways: the program Owens.f was used by two 
members of our group to perform independent profile fits, and a third member of our 
group performed both a profile fit using the code of Fitzpatrick (Spitzer & Fitzpatrick 
1993) and a curve-of-growth analysis. Each analysis was performed with different choices 
of continuum regions to fit, different choices of polynomial order for the continuum fits, 
and different combinations of background level and LSF. The resulting column densities, 
with 2 a uncertainties, were: log Noi = 14. log Noi = 14.48+°$, log N 0 i = 
14.50 ± 0.04, and log N OI = 14.5ll{j;o6. The corresponding b-values and 2 a uncertainties 
from each fit were: 2.8 ± 0.2, 3, 2.5 ± 0.4, and 2.7±0.3. The b-value for the second fit 
is poorly-determined because only unsaturated lines were included, which are insensitive 
to b. The first two b-values were from simultaneous fits to O I, N I, and N II lines in the 
FUSE data; if only the O I lines are included in the first fit the uncertainty increases to 
±0.40 km s -1 . The b-values fall somewhat below the value predicted from the line of sight 
velocity structure derived from the GHRS data, but are consistent within 2 a. The above 
determinations are in good agreement with one another; we adopt the value log N 0 i = 
14.49 ± 0.08 (2 a). A sample of the fits to the O I lines is shown in Figure 6, and the curve 
of growth is shown in Figure 7. 

A number of tests were performed to examine the magnitude of the various possible 
systematic effects described in Section 4.2. Increasing the background level by a factor of 
2 in the SiC channels and a factor of 4 in the LiF channels increased the value of log N 0 / 
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by 0.02 dex, decreasing the background by a factor of 2 decreased log Noi by < 0.01 dex. 
Fitting only detector 1 data or only detector 2 data, to test for fixed-pattern noise effects 
and detector-dependent r solution or background effects, changed log Noi by < 0.01 dex. 
Eliminating any single line from the fit, to test for erroneous f- values or fixed-pattern noise 
effects, had no effect on log Noi ■ Eliminating both the line at 1039A and the multiplet at 
988A from the first profif 1 fit causes log Noi to increase by 0.04 dex, but the uncertainty 
of the result doubles. Eliminating these lines from the COG analysis causes the resulting 
log Noi to increase by only 0.02 dex. The LSF was varied by doubling the amplitude of 
the broad component, which increased log Noi by 0.02 dex, and eliminating the broad 
component (with a corresponding increase in the background level to match the cores of 
the saturated Lyman lines), which decreased log Noi by 0.02 dex. These effects have all 
been included in the uncertainties quoted above for the first profile fit. 


4-3.4 ■ N I Column Density 

As with O I, the FUSE bandpass contains numerous Nl absorption lines of varying 
strength. In this case the low column density to HZ43A means that many of the lines that 
are ordinarily the most useful are not detected. High quality data were obtained in all 
channels in the MDRS observation, and usable data were obtained in all channels in the 
LWRS observation. We analyzed both the complete dataset and the pure night data for the 
LWRS observation and found that the geocoronal emission from Nl during the day portion 
of the orbit did not significantly affect the absorption line profiles. The lines strong enough 
to be useful in the FUSE data were the triplet at 1134A and the lines at 953.655A and 
953.970A. 

The MDRS data were analyzed by having two members of the group perform 
independent profile fits and a third member perform a curve-of-growth fit. The curve-of- 
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growth analysis also included the N I multiplet 1199. 550A, 1200. 223A, and 1200. 710 A from 
the GHRS G160M dataset. For the MDRS data, the results were log N n1 = 13.51 1q'o 6> 
log N ni = 13.51^q q 5, and log N ni = 13.51^ 05 (2 cr uncertainties). The b-value and 2 a 
uncertainty derived from the curve of growth analysis was S.O+J^km s -1 ; no determination 
of b was possible from the profile fits because the lines used were all on the linear part of 
the curve of growth. The adopted value for log N^i is 13.51 ± 0.06. The fit to the MDRS 
N I profiles is shown in Figure 8, and the curve of growth is shown in Figure 9. 

The N I fits were subjected to the same tests for systematic errors as described in the 
previous section for O I. Because all of the N I lines in the FUSE band relatively weak, they 
are less sensitive to variations in the LSF or the background than partially saturated lines; 
none of these tests caused a variation in the log N ni greater than 0.01 dex. 


4-3.5. D I Column Density 

Because of the low H I column density along this line of sight, absorption by D I is 
expected to be detectable at Lyman-a, Lyman-/?, and possibly at Lyman-y. The Dl 
Lyman-a profile is shown in Figures 2 and 4; it is clearly well-resolved. The D I Lyman-/? 
absorption is clearly seen in all FUSE channels in the MDRS observation, but is clearly 
detected only in the night portion of the LWRS observation in the LiFl channel. The 
MDRS SiC2 Lyman-/? spectrum was rather noisy, however, and was not included in the fits. 
Dl Lyman-7 absorption is not detected in either SiC channel. 

The GHRS Lyman-a profile was first fit independent of the FUSE data. The resulting 
value for the column density was log Npi = 13.16 ± 0.04 (2 cr), with a b-value of 
7.1 (implying the line is well-resolved by the GHRS echelle resolution). Changing the 
background level from our best estimate by ± 2% of the continuum level changed log N& j 
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by ±0.02 dex. A simultai eous fit to the GHRS Lyman-c* spectrum and the FUSE MDRS 
Lyman-/? data from channels LiFl, LiF2, and SiCl gave log Npi = 13.15 I 0.045 a )- 

Varying the fitting parameters for the FUSE data as described in Section 4.3.3 above for 
Oi resulted in variations in log Npj of ±0.01. Profile fits to the FUSE MDRS data alone 
gave log N di = 13.11 ± 0.20. The LWRS LiFl Lyman-/? spectrum was not included in the 
fits. The D I line is detected in this spectrum, and is consistent with the fits to the other 
data, but the uncertain corrections for airglow contamination in the LWRS spectrum, even 
for the night only data, means that these data will not significantly constrain the overall fit. 
The FUSE D I line profiles are shown in Figure 10, along with the best fit to the combined 
FUSE and GHRS datasets. 


4-3.6. Other Species 


A number of other species were detected in the ISM along the line of sight to HZ43A. 
These include C II, Cm, Nil, Ovi, Sin, Ari, and possibly Nm in addition to the Mgii 
and Fell mentioned earlier. The column densities for these species are listed in Table 5. 

Of particular interest is t lie column density for Nil (log iVv//= 13. 62 +^ 7 ), which slightly 
exceeds that of N 1 . The ionization potential for N I is greater than that for H I, so ordinarily 
one expects N to be predominantly neutral. The local ISM is dominated by hot gas that 
may be far from ionization equilibrium, however, so estimates of the total abundance of 
nitrogen in the local ISM must account for possible large ionization corrections to the 
measured N I abundances. The N II transitions at 915.613A and 1083. 993A are ordinarily 
strongly saturated, but the total column density to HZ43A is low enough that an accurate 
N II column density can be determined. The N II absorption lines and best-fit line profiles 
are shown in Figure 11. The N ill column density quoted in Table 5 is an upper limit 
because the N ill 989.799A is blended with the Sill 989.873A line. If the Sin column density 
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derived from the 1193. 290A line and the b-value derived from Ol are used to predict the 
equivalent width of the Si II 989.873A line, then it is likely that at least half, and perhaps 
all, of the measured equivalent width at 989. 8A is caused by Sill rather than Nlll. Even if 
the measured absorption results primarily from N III, the upper limit suggests that N ill is 
a minor contributor to the total N column along this sightline. 

No H 2 was detected. The range of possible column densities was determined by using 
the profile-fitting program Owens to examine the regions of the spectrum in the vicinity of 
the two-dozen strongest H 2 absorption lines that were not blended with any other features. 
The column density of H 2 was then gradually increased until the change in chi-squared was 
unacceptable. The resulting 3 — a upper limit to the column density of H 2 was 1.2 x 10 13 
cm -2 . 

The column density for C II shown in Table 5 is larger than that for O I, possibly 
suggesting that C is over-abundant on this sightline, that O is depleted, or that there is 
substantial component of H II gas along this line of sight. The large column density of 
N II supports the latter explanation. It should be noted also that the column density for 
C II is particularly uncertain, because the saturation effects are large and depend on the 
assumption that the b-value for C II is the same as that of O I. If there is a substantial 
contribution to the column density from H II gas, the effective b-value for C II will likely 
be larger than that for O I and the C II column density will be overestimated accordingly. 
This source of systematic error is not included in the uncertainties quoted in Table 5; the 
column density derived from the apparent optical depth method (log N(Cll) = 13.63) is a 
reasonable lower limit. 

We note that Cll* at 1037.018 A is not detected in our FUSE observations. The 30- 
upper limit to the equivalent width of this line in the LiFlA data is W\ < 5.5 mA, including 
a crude estimate of the continuum placement uncertainties. This implies log N ( C iT) < 12.7 
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(3cr) along the HZ43A sight line. If one assumes the excitation of Cn* is due to electrons 
the electron density may be calculated by comparing the column densities of the excited- 
and ground-state C II: 


n e 


0.18 ^“ V ' 2 

JV(Cii) 


(see Spitzer & Fitzpatrick 1993). Assuming T = 5353 K (see Section 4.3.1) and that 
log N(C il) = 14.83 (derived using the Ol 5-value), we find n e < 0.1 cm -3 . Given the 
uncertainties in the latter assumption, however, a more conservative value is n e <C 1.5 cm 
derived using the lower limit for the Cn column density determined with the apparent 
optical depth method. 


5. Summary 

The line of sight to HZ 43A was investigated using a combination of FUSE spectra, 
archival GHRS spectra, arid results obtained from EUVE spectra taken from the literature. 
The line of sight to HZ43A is particularly simple, with only a single velocity component 
discernible. An additional high-temperature component was found with a very low column 
density, detectable only in Hi. The line of sight to HZ43A exits the LIC after a very short 
distance, and the velocity of the main component of the absorbing gas is not compatible 
with the velocity of the nearby G cloud in the LISM, hence this material is presumed to 
be located in the North ( lalactic Pole cloud. Doppler widths for each species determined 
from profile fits or curve of growth analyses of the FUSE data were consistent with the 
temperature and turbulent velocity for the gas derived from the GHRS high resolution data. 
Consistent fits to the Dl column density were obtained when fitting the FUSE and GHRS 
data separately and simultaneously, indicating that the D I Lyman-c* line is not affected by 
saturation effects. The F USE spectra were used to determine O I and D I column densities 
by both profile fitting and curve of growth analyses, with consistent results obtained in each 
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case. Careful attention was paid to numerous potential sources of systematic error. Our 
adopted Hi column density, log N H1 = 17.93 ± 0.06 {2 a), is the mean of our measurement 
from the GHRS Lyman-a profile and the EUV measurements of Dupuis et al. (1995), 
Barstow et al. (1997), and Wolff et al. (1999). 

Our results are summarized in Table 6. All uncertainties in the table are 2 a (but 
values quoted below from the literature are 1 a). Our result for D/H along the line of sight 
to HZ43A is (1.66 ± 0.28) x 10 -5 . The GHRS dataset analyzed in this work has been 
previously analyzed by Landsman et al. (1996), who obtained the similar value for D/H of 
1.6 x 10 -5 . Our value for D/H along this sightline exceeds by approximately 1.5 a both the 
mean value reported by Linsky (1998) for measurements of D/H within 100 pc of the sun, 
(1.47 ± 0.10) x 10 -5 , and the mean value of the FUSE measurements reported in this first 
set of papers (Moos et al. 2001). 

Our result for Ol/H I is (3.63 ± 0.84) x 10 -4 . Converting this ratio to a logarithmic 
abundance log (Ol/Hi) + 12.00 gives 8.56+oJj. The ionization balance for both O and H 
are linked by resonant charge exchange reactions (see Jenkins et al. 2000 and references 
therein), and there is no apparent H2 present, so Ol/Hl should be representative of O/H 
along this line of sight. Meyer et al. (1998) report an average gas-phase abundance for 13 
sight lines of (3.43 ± 0.15) x 10 -4 (after correction for our preferred Ol A1355 oscillator 
strength of / = 1.16 x 10~ 6 from Welty et al. 1999), and they estimate no more than 
0/H~ 1.8 x 10 -4 can be incorporated into grains, implying an upper limit to the total 
(gas+dust) ISM oxygen abundance of 5.2 x 10 -4 . Our O/H measurement is consistent with 
the average gas-phase value of Meyer et al. (1998). 

Sofia Sz Meyer (2001b, a) have demonstrated that the value of O/H determined from F 
and G star photospheres, (4.45 ± 1.56) x 10 -4 , is consistent with two recent determinations 
of the solar abundance, (5.45 ± 1.00) x 10"" 4 (Holweger 2001) and (4.90+053) x 10~ 4 (Allende 
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Prieto et al. 2001), and with the total oxygen abundance estimated from the work of 
Meyer et al. (1998). Sofia k Meyer (2001b, a) argue that the recent solar system values 
should therefore be used is the ISM abundance standard for oxygen. Our value for O/H 
in the gas phase along the line of sight to HZ 43 A is significantly lower than each of these 
values, although not by more than 2 a (given the large uncertainties in the solar system 
determinations). Our result is in agreement with the compilation of B-star abundances 
presented by Sofia k Meyer (2001b, a), (3.50 =t 1.33) x 10“ 4 , although these authors argue 
that B-star abundances are unlikely to be appropriate measures of the ISM standard. It 
seems likely that the gas- phase oxygen abundance measured by Meyer et al. (1998) implies 
modest depletion of oxygen in the ISM towards more distant stars; given the agreement of 
our O/H measurement with the Meyer et al. average, it seems likely a similar degree of 
oxygen depletion is present in the cloud observed along the HZ43A sight line. 

The ionization fractions for N and H are not as strongly coupled as those of O and H 
(Jenkins et al. 2000), so i’ is possible for N to be more highly ionized than H, and therefore 
for the gas-phase abundance of N relative to its solar system value to be less than that of 
O. If we combine our column densities for N I, N II, and the maximum likely column density 
for N ill (assuming only half of the measured absorption is due to Sin), we obtain a total 
nitrogen column density of = 7.88t};43 x 10 13 > an( i a ratio N/Hl = (9.26 ± 1.87) x 10 5 . 
This ratio provides an upper limit to N/H, corresponding to the limit that that hydrogen 
along this line of sight is predominantly neutral. Given the significant column density of 
N II, it is reasonable to expect that there will also be a significant column density of H II 
present. Converting to a logarithmic scale with log N H j = 12.00 gives an abundance for N 
of 7.97^0 i 0 ’ This is consistent with the solar abundance value of Grevesse k Noels (1993), 
7.97, and the more recent value from Holweger (2001) of 7.93 ± 0.11, and somewhat exceeds 
the value of 7.81 for nearby B stars (Gies k Lambert 1992). We can set a lower limit to 
the N abundance by setting the ionization fraction of H to the upper limit of 50% found 



-40 - 


by Dupuis et al. (1995) (the other EUVE analyses obtained lower ionization fractions for 
hydrogen). The resulting logarithmic abundance is 7.79louo> essentially equal to the value 
for B stars, and consistent with the ISM gas phase abundance of Meyer et al. (1997). Thus 
N is not significantly depleted onto dust grains along this line of sight. 

Further discussion of the implications of these results and comparisons with other 
sightlines can be found in the companion paper by Moos et al. (2001). 

This work is based on data obtained for the Guaranteed Time Team by the NASA- 
CNES-CSA FUSE mission operated by the Johns Hopkins University. Financial support to 
U. S. participants has been provided by NASA contract NAS5-32985. French participants 
are supported by CNES. We would like to thank Ed Fitzpatrick for use of his profile fitting 


software. 
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Relative Velocity (km/s) 


Fig. 1. — Absorption profiles of Dl, Fell, and Mgll are plotted as a function of heliocentric 
velocity. The histograms are the GHRS echelle measurements, and the smooth line is the 
single-component fit. No evidence for a second velocity component is seen. The Dl line sits 
on the broad wing of the interstellar H I Lyman-a profile; the truncated peak at +65 km s 1 
in the D I plot is geocoronal Lyman-a emission. 
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Wavelength (A) 

F ig 2 — The GHRS eclielle spectrum of the Lyman-a profile is shown as the histogram, 
the stellar model as the smooth curve. The model is shifted to a heliocentric velocity of 
+21 km s -1 . The positions of the interstellar Hi and Dl absorption lines are marked. The 
truncated peak at 1215.6.4 marked with an earth symbol is geocoronal H I Lyman-a emission. 
The GHRS spectrum was normalized to match the mean flux of the model over the regions 
1212.25A-1213.25A and 1218.0A-1218.5A; additional low-order polynomial corrections were 
applied to the continuum flux when fitting (see text). The zero level is shown as the dashed 


line. 
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Fig. 3. — Plots of the x 2 contours are shown for fits to the GHRS Lyman-a profile as a 
function of log TV h / and photospheric velocity. The top plot shows the results for 4 th order 
polynomial continuum fits, and the bottom plot shows the results for 7 th order polynomial 
continuum fits. 
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Fig. 4 — The normalized GHRS spectrum of the Lyman-a profile is shown as the histogram, 
and the best fit is shown as the smooth curve. The thin curves show the continuum fit and 
the individual H I and D [ profiles; the dashed curve shows the low-column density high- 
temperature Hi component; and the thick solid line shows the overall fit. The stellar model 
in this fit was shifted to a lieliocentric velocity of 26 km s -1 , which provided slightly better fits 
than the nominal value of 21 km s -1 . The positions of the interstellar H I and D I absorption 
lines are marked. The truncated peak at 1215. 6A marked with an earth symbol is geocoronal 
Hi Lyman-a emission. The zero level is shown as a dashed line. 




52 - 



0 I AUijxzntlZ llKMdjBSSai 1 ■ i I I 1 1 1 1 ■ 1 1 1 1 ■ I — 

912 913 914 915 916 


Wavelength (A) 

Fig. 5. — The FUSE MDRS SiClb spectrum of the converging H I Lyman series is plotted 
as the histogrammed points, the profiles resulting from the best fit to the Lyman-a profile 
as the smooth line. The assumed continuum is shown as the thin solid line. The LWRS 
spectrum is essentially identical, but with lower signal to noise. 
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Fig. 6. — The observed FUSE MDRS O I line profiles are shown as the histogrammed data 
points, the continuum fits are the thin solid lines, and the profiles computed from the adopted 
Ol column density and doppler width is plotted as the thick solid curve. The data are un- 




1039 1039.5 


binned. 
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Fig. 7. — The equivalent widths of Ol lines measured from FUSE MDRS data are plotted, 
along with the best-fit turve of growth (b = 2.72km s" 1 ). The inset shows contours of 
constant A X 2 in the b~lcgN plane. The measurements of lines from multiple channels have 
been offset slightly in log NAf to avoid confusion. The wavelengths of the lines are marked 
along the top of the figure. The plotted errorbars indicate the la measurement uncertainties 

in the equivalent widths. 
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Fig. 8. — The FUSE MDRS N I line profiles are shown as the thin histogram and the best 
fit is shown as the thick solid line. The continuum fit is the thin smooth line. The SiC data 
have been binned by two pixels for plotting purposes. 
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Fig. 10. — The FUSE spectra of Hi and Di Lyman-/? and Lyman-7 absorption are plotted 
as the thin histogram, and the best fit to the combined GHRS and FUSE data as the thick 
smooth line. No D I absorption is seen at Lyman-7, but it is clearly seen at Lyman-/?. The 
geocoronal Lyman-/? emission is clearly seen in the red wing of the interstellar absorption; it 
is well-separated from D 1. Geocoronal emission at Lyman-7 is present, but quite weak. The 
O I 1026. 473A absorption predicted from the measured O I column density is also plotted; 
no absorption is seen at this wavelength so it has been assumed that the oscillator strength 
for this line is in error. It was not included in any fits. 
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Pig 11 The FUSE MDRS N n line profiles are shown as the thin histogram and the best fit 

is shown as the thick solid line. The continuum fit is the thin smooth line. The neighboring 
interstellar Hi line profiles are labelled with tick marks in the plot of the 915.663Aline. 
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Table 1. Target Summary for HZ43A 


Quantity 

Value 

Reference 

Spectral Type 

DAI 

i 

c* 2 ooo a 

13'*16 m 21 s 78 

2 

^2000 a 

+29°05'55"5 

2 

l 

54.11° 

2 

b 

+84.16° 

2 

d h (pc) 

68 + 13 

3 

V 

12.914 

2 

B-V 

-0.31 

2 

T eff (K) 

50,900 

1 

logy (cm s -2 ) 

8.0 

1 


a Corrected for proper motion to epoch 
2000.5 

b Trigonometric parallax 


References. — (1) Finley et al. 1997; (2) 
Bohlin et al. 1995; (3)van Altena et al. 1995. 
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Table 2. Log of FUSE observations. 


Dataset ID 

Aperture 

Obs date 

Texp (total)® 

Texp 

N cxp h 

M1010501 

LWRS 

Feb 19 2000 

6092 

2333 

8 

P1042301 

LWRS 

Apr 22 2000 

14447 

5136 

23 

P1042302 

MDRS 

Feb 08 2001 

39606 

20572 

85 


a Exposure luration in seconds. 


b Number of individual exposures. 


Table 3. 


Log of GHRS observations. 


Dataset 

ID 

Grating 
and order 

Aperture and 
substep pattern 

Wavelength range 

A 

Start 

Time (UT) a 

Exposure 
Time (s) 

Z2R50104T 

Ech A - 46 

LSA 7 

1212-1218 

00:21 

7834 

Z2R50107T 

Ech A - 46 

LSA 7 

1212-1218 

05:05 

7181 

Z2R50204T 

Ech B - 24 

LSA 7 

2376-2387 

09:59 

2176 

Z2R50205T 

Ech B - 20 

LSA 7 

2791-2807 

11:27 

2176 

Z2R50209T 

G160M 

SSA 5 

1192-1228 

13:10 

1414 


a All exposures began on July 30 1995. 
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Table 4. Equivalent Width Measurements 


Species A c a logA/ b Channel W x [mA] c S/N d 


FUSE 


CII 

C III 977.020 2.870 

N I 953.655 1.376 

N I 953.970 1.521 

N I 1134.165 1.238 

N I 1134.415 1.528 

N I 1134.980 1.693 


LiFlA 

37.81J;? 

42 

LiF2B 

39-8 tU 

26 

SiClA 

39.4+2 4 

19 

SiClB 

14.6+2:2 

16 

SiC2A 

15-611:1 

22 

SiClB 


15 

SiC2A 

4-7 + _li 

17 

SiClB 

9-6l^ 

16 

SiC2A 

n 

' *^-2.6 

20 

LiFlB 

r: c+1.7 

47 

LiF2A 

5-4+Jg 

47 

LiFlB 

8.9^1; 7 0 

46 

LiF2A 

9.2l+| 

43 

LiFlB 

14.1^; 6 5 

46 

LiF2A 

iiilui 

43 


1036.337 2.088 



Table 4 — Continued 


Species 

A c a 

log A/ b 

Channel 

W x [mA] c 

S/N d 

N II 

915.613 

2.180 

SiClB 

20.3 tfl 

13 

N II 

1083.994 

2.097 

SiC2B 

23.3l|;? 

14 

N III 

989.799 

2.085 

SiC2A 

< 8.2l?;? e 

25 

01 

936.630 

0.534 

SiClB 

o q-t- 2.9 
0.0_2 7 

13 




SiC2A 

11-411;? 

15 

01 

948.685 

0.778 

SiClB 

7-811;? 

12 




SiC2A 


15 

01 

971.738 

1.123 f 

SiClB 

I5.9ll;g 

10 




SiC2A 

16-5111 

13 

01 

976.448 

0.509 

SiClB 

8-411:4 

15 




SiC2A 

8011 q 

18 

01 

988.655 

0.914 

SiClB 

16.011:1 

20 




SiC2A 

i3oll;I 

23 




LiF2B 

15-411;? 

14 

01 

988.773 

1.662 

SiClB 

27.51H 

20 




SiC2A 

30.411;? 

24 




LiF2B 

26.il?;l 

14 
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Table 4 — Continued 


Species 

A c a 

log A / b 

Channel 

W x [mA] c 

S/N d 

01 

1039.230 

0.980 

LiFlA 

19. 11^2 

42 




LiF2B 

16.9l^ 

26 




SiClA 

16.4±i;| 

19 

0 VI 

1031.926 

2.137 

LiFlA 

8 - 5 ^ 1 .7 

39 




LiF2B 

8-411:1 

25 

Ar I 

1048.220 

2.442 

LiFlA 

o 0 + 1-2 
^-^- 0.8 

49 




LiF2B 

< 2.4 (2a) 

25 

GHRS 

N I 

1199.550 

2.192 

G160M 

27.6111 

17 

N I 

1200.223 

2.015 

G160M 

2 O. 2 I 44 

17 

N I 

1200.710 

1.713 

G160M 

n.olil 

16 

Mg II 

2796.352 

3.236 

Ech-A 

55l 3 3 

12 

Mg II 

2803.531 

2.933 

Ech-A 

37.911:2 

13 

Si II 

1193.290 

2.775 

G160M 

24.711S 

16 
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Table 4 — Continued 


Species 

V 

log A/ b 

Channel 

W x [mA] c 

S/N d 

Fe II 

2382.765 

2.882 

Ech-A 

20.511;* 

12 


a Central wavelength of each transition from Morton 
2001 . 

b Adopted product of wavelength and oscillator strength 

for each transition, taken from Morton 2001. 

c Measured equivalent width and la uncertainties (in 

mA). The uncertainties contain contributions from both 

statistical and systematic error sources. 

d Empirically estimated signal-to-noise ratios (per un- 
binned pixel) for continuum regions near the absorption 
lines of interest. 

e The N III absorption contains a substantial contribu- 
tion from Si II at 989.873 (« +21 km s -1 from the center 
of the N III line). The equivalent width quoted may be 
dominated by the latter, but represents a firm upper limit 
for the amount of N III along the HZ 43A sight line. 

f This transition is a blend of three O I lines, the sum of 
whose oscillator strengths are given here. 
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Table 5. Curve of Growth Column Densities 


Species 

log AT ± 2 cr a 

b±2a h 

A log N c 

C 11 

14-83 t° Q £ 

d e 

+1.20 

C HI 

12.59t° 0 ;° 9 0 

d 

+0.15 

N I 

13.5llo.o5 

3.01J:? 


N II 

13.6212:1? 

d 

+0.32 

N III 

< 12.981 ° 0 % 

d 

+0.07 

0 I 

1 4 M + 007 

0 70^“^ -33 
* z -0.26 


0 VI 

i2.9HS:SS 

oo f 


Si II 

12.8512:1? 

d 

+0.21 

Ar 1 

ii-96l2i2 

d 

+0.00 

Mg II 

12.41±0.02 

2.8±0.3 


Fe II 

12.17±0.04 

2.78±1.15 



l Log of column density derived from a curve 
of growth fit, with estimated 2cr uncertainties. 
For most species not enough transitions were 
observed to independently fit a curve of growth 
and determine a b- value. In these cases we have 
adopted the b - value of 2.72l2?2 (2cr) from O I. 

b Doppler parameter in km s 1 with esti- 
mated 2 a uncertainties. 

°For those species for which we have adopted 

the b - value from O I we give the amount by 
which the curve-of-growth-corrected column 
density given in column 2 of this table exceeds 
the average integrated apparent column den- 
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Table 6. HZ 43 Column Densities and Ratios 11 



Hi 

Dl 

Oi 

Nl 

Log N 

17.93 ± 0.06 

10 1 r +0.040 
lO.lO _ 0 045 

14.49 ± 0.08 

13.51 ± 0.06 

Ratios 

Di/Hi 

Ol/Hl 

Ni/Hi 

Di/Oi 

Dl/Nl 

1.66 ±0.28 x 10“ 5 

3.63 ±0.84 x 10" 4 

3.80 ± 0.74 x 10~ 5 

4.57 ±0.96 x KT 2 

4.37 ±0.74 x 10" 1 


a all uncertainties are 2 o 



